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It is shown that when the nanobubble contact line is pinned to a penetrating tip the interface 
behaves like a Hookean spring with spring constant proportional to the nanobubble surface tension. 
Atomic force microscope (AFM) data for several nanobubbles and solutions are analysed and yield 
surface tensions in the range 0.04-0.05 N/m (compared to 0.072 N/m for saturated water), and 
supersaturation ratios in the range 2-5. These are the first direct measurements of the surface tension 
of a supersaturated air-water interface. The results are consistent with recent theories of nanobubble 
size and stability, and with computer simulations of the surface tension of a supersaturated solution. 


I. INTRODUCTION 

Experimental evidence for nanobubbles was first pub¬ 
lished in 1994,— and since then their existence has 
been confirmed from various features of the mea¬ 
sured forces between hydrophobic surfaces includ¬ 
ing a reduced attraction in de-aerated water^^— and 
from images obtained with tapping mode atomic force 
microscopy— For a recent review of theory and 
experiment, see Ref. M 

The initial controversies over the existence of nanobub¬ 
bles —that they should have an internal gas pressure of 
10-100 atmospheres that would cause them to dissolve in 
microseconds— have largely been resolved. First it was 
shown that nanobubbles can only be in equilibrium in 
water supersaturated with air^iS and then it was shown 
that the surface tension of a supersaturated solution must 
be less than that of a saturated solution;^ Finally, it 
was shown that nanobubbles with a pinned contact rim 
are thermodynamically stabled This means that con¬ 
tact line pinning is a necessary and sufficient condition 
for nanobubbles to be in mechanical and diffusive equi¬ 
librium with the supersaturated solution. 

The significance of the reduction in surface tension is 
that the internal gas pressure, which can be calculated 
from the Laplace-Young equation, is much less than those 
initial estimates used to argue against nanobubbles. It 
also means that the degree of supersaturation of the solu¬ 
tion, necessary for diffusive equilibrium of the nanobub¬ 
ble, is reduced to realistic levels that are attainable in 
the fluid cell. 

For many it is surprising that purely on thermody¬ 
namic grounds (ie. no additives or surfactant) the surface 
tension for nanobubbles should be reduced from the usual 
value of the air-water interface. Nevertheless, this result 
is firmly established by thermodynamics^ density func¬ 
tional theory,— and computer simulation.The 
result is not widely acknowledged within the nanobub¬ 
ble field, possibly because many practitioners place more 
trust in experimental measurement than they do in ther¬ 
modynamics or in mathematical equations. To close this 
gap, it would be desirable to measure directly the surface 
tension of nanobubbles. 

Measuring the surface tension of the supersaturated 



FIG. 1: Hemispherical bubble penetrated and deformed by a 
conical tip (blunt radius ro, height zq) upon which it sticks 
at r\. 


air-water interface is a worthwhile experimental goal that 
has application beyond nanobubbles. For example, in at¬ 
mospheric physics the nucleation of cloud droplets always 
occurs from a supersaturated atmosphere, and so the rate 
of change of surface tension with the degree of supersatu¬ 
ration is a key input necessary for quantitative modeling 
in that field. 

A supersaturated solution in the current context means 
that the concentration of air in the water is higher than 
it would be if the water were equilibrated with air at the 
current temperature and pressure. Supersaturation can 
be readily achieved by previously equilibrating the water 
with air at a higher pressure or at a lower temperature. 
The difficulty in measuring directly the surface tension 
of a supersaturated solution is that the change in surface 
tension is determined by the concentration of air within 
the first nanometer of the interface, but, due to diffu¬ 
sion across the interface, this region is always saturated 
rather than supersaturated, at least if the measurement 
is performed on a macroscopic droplet or bubble exposed 
to the atmosphere. In the case of nanobubbles, how¬ 
ever, the solution is supersaturated immediately in the 
vicinity of the interface due to the high internal gas pres¬ 
sure of the nanobubble, and so the surface tension of a 
nanobubble must be that of a supersaturated solution. 
If one could measure the surface tension of the nanobub¬ 
ble, then one would have a means to determine the rate 
of change of surface tension with supersaturation. The 
present author knows of no previous measurements of the 
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surface tension of the supersaturated air-water interface. 

The present paper is concerned with modeling 
atomic force microscope (AFM) force measurements on 
nanobubbles when the tip of the cantilever penetrates the 
nanobubble. The aim is to relate the surface tension of 
the nanobubble to the measured force. The results are 
applied to force data from several nanobubbles. In all 
cases analyzed, the surface tension was always less than 
that of the saturated air-water interface, and the super¬ 
saturation ratio of the solution was always greater than 
unity. 


II. PENETRATED HEMISPHERICAL BUBBLE 

Figure[I]is a sketch of a hemispherical bubble on a solid 
surface that is deformed by a penetrating conical tip. 
In a cylindrical coordinate system, the bubble contacts 
the surface at (r3,0) and the tip at {ri,zi). The bubble 
is pinned or fixed at these two radii. The blunt tip of 
the cone is at a height (separation) above the substrate 
of zq for r < vq. The half angle is given by tana = 
{ri — ro)/(zi — Zq). The case of a positive load, F > 0, 
is shown in the figure; for a negative force the bubble is 
extended and the dimple vanishes. 

The following analysis is similar to that given earlier 
in that it is based on a small force expansioni^ The ear¬ 
lier analysis explicitly included the effects of a surface 
or interaction force between the probe and the bubble, 
whereas here the bubble is penetrated by the probe and 
makes either stick or slip contact with its sides. The ear¬ 
lier analysis was for a bubble mobile on the substrate and 
for fixed number of air molecules, whereas here the bub¬ 
ble is pinned at rs and is in diffusive equilibrium with 
the supersaturated solution. Due to the present pinned 
contact rim, the solid surface energies and the contact 
angle play no role in the present analysis, whereas they 
did for the case of a mobile bubble analyzed in Ref. 1^. 
These differences in the model make the present analysis 
considerably simpler and shorter than that of the earlier 
case. Despite these differences in the model, the qualita¬ 
tive conclusion is the same in both cases, namely that the 
bubble interface behaves as a Hookean spring. The quan¬ 
titative expression for the spring constant found here is of 
course specific for the present model of a bubble pinned 
at rs. 


A. Undeformed Bubble 


As recently shown^ once the bubble is pinned at the 
contact radius ra, it is thermodynamically stable at the 
critical radius and density. Accordingly, the undeformed 
bubble is the critical bubble. It is of course hemispheri¬ 
cal, and its radius is the critical radius^ii^ 


27 27 

(s - l)p Ap ■ 


( 2 . 1 ) 


Here s > 1 is the supersaturation ratio, (the solution is 
necessarily supersaturated with air in order for the bub¬ 
ble to be in diffusive equilibrium) ,— 7 is the liquid-vapor 
surface tension (of the supersaturated interface) p 

is the pressure of the reservoir (taken to be atmospheric), 
Ap is the excess pressure of the bubble, fee is Boltzmann’s 
constant, and T is the temperature. For simplicity, the 
reservoir pressure and the saturation vapor pressure are 
approximated as equali^ 

The undeformed bubble has volume^! 


14 = TT [Rc^c - ^c/3] , 
liquid-vapor surface area 


Ac = 2 'kRc 
and apex height 


Rc - ^/i?2 - ri 


= 27ri?c^c, 


Zc = Rc - \/Rc - rl- 


The undeformed profile is 

Zc{r) = Zc- Rc + 


( 2 . 2 ) 


(2.3) 


(2.4) 


(2.5) 


The undeformed bubble contains Nc = spYc/k^T gas 
molecules. In the present case of fixed contact rim r 3 , 
thermodynamic stability holds simultaneously for num¬ 
ber and volume fluctuationsHence one does not have 
to insist upon constant number as the author has had to 
do in previous work.i^ii^ 


B. Deformed Bubble: Prick Stick 


The bubble profile is z{r), which ends at the contact 
points (ri,zi) and (r3,0). The contact points on the 
blunt tip satisfy ri — tq = [zi — zq] tana, where a is the 
cone half angle, and zq is the height of the tip of the tip 
above the substrate, which is also called the separation. 
The tip is taken to be perfectly blunt, which is to say that 
its end is a disc of radius tq at zq- Although in principle 
one could also have contact at [ri,ZQ) for ri < rp, this 
case will be excluded in the numerical results below. A 
perfectly sharp tip has rg = 0. 

The volume is 


/•rs 

V[z] = 2 tt drrz{r)+TTrlzi — Vt 

Jri 

1 27rrf 

= 27r / arrz(r) + 

J n 


3 tan a 


■ Tirl 


zo - 


ro 


tana 


' 0 


3 tan a 


The volume of the blunt tip inside the bubble is 


Vt = 


TrrV 


3 tan a 3 tan a 


( 2 . 6 ) 


F, 


(2.7) 
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The terms irr^zi — Vt give the vapor volume beneath 
the tip, r < ri. These may be neglected for the profile 
differentiation since they are constant. The fluid-vapor 
interfacial area is 


rr3 

A[z] = 27r / dr ri/l -|- z'(r)^ 
J ri 


( 2 . 8 ) 


Also 


SA[z] = 


dr 1 — — ^7'{r) 




n 1 d 

—27r / dr 
Jri 


V^ + z'{r)^ 

z'ir) 


dr 


y/1 + z'ij)^ 


(2.15) 

i5z(r), 


The total entropy of the bubble and reservoir is a func¬ 
tional of the profile and a function of the other thermo¬ 
dynamic parameters, 


^totN = 5b(iV,F,r)-lA[z]-|l-[z] + |fV(2.9) 

The solid-air and solid-liquid surface energies (for both 
the substrate and the tip) are constant because of the 
fixed contact radii and so are neglected here. Here 
/j, = k-QT \n[spK^/k-QT] is the chemical potential of the 
air. This expression is appropriate for the case that the 
bubble can exchange number, volume, and area with the 
reservoir. 

It is assumed that thermal equilibrium holds. The en¬ 
tropy of the bubble may be taken to be that of an ideal 
gas for thermal equilibrium. 


5b(iV,H,r) 


k^N 


1-ln 


NM^' 

V 


( 2 . 10 ) 


where A is the thermal wave length. One could in 
fact use the entropy of a real gas instead of this since 
the only things that enter below are its thermody¬ 
namic derivatives, pb = TdS\y{N,V,T)/dV, and pb = 
-TdSh{N,V,T)/dN. 

Obviously, setting the number derivative of the total 
entropy to zero yields the equilibrium condition Pb = 
which for an ideal gas is the same as 


following an integration by parts and the vanishing of the 
perturbation at the boundaries. From this one has 


5A[z\ ^ d rz'{r) 

5z{r) ^ dr + z'{rY 


(2.16) 


Inserting these into the functional derivative of the total 
entropy and setting the latter to zero gives a differential 
equation (the Eular-Lagrange equation) for the optimum 
profile. 


0 = 27rr Ap -|- — 
dr 


2TTrz'(r) 
V'l -bz'(r)2 


(2.17) 


Due to diffusive equilibrium, the excess pressure is a con¬ 
stant, Ap = (s — l)p. 

The first integral of this is 

const. = - ( 2 . 18 ) 

" v'TTvyy' 


1. Force 

Including a spring attached to the cantilever with 
spring constant fct, the extended total entropy is 

Stot,k = 5'tot([-z],ri) - - Zcf. (2.19) 


( 2 - 11 ) 

Invoking the usual variational properties of equilibrium 
thermodynamics^ the number is fixed at this value in 
all that follows. 

The functional derivative of the total entropy with re¬ 
spect to the bubble profile is 

(5gtot[^] _ Ap (5y[z] 'y 5A[z\ 

5z{r) T 6z{r) T 5z{r) 

where Ap = p\^ — p = (s — l)p is the excess pressure of 
the bubble. One has 

rrs 

SV[z] = 27r / drrSz(r), (2.13) 

J n 


from which it follows that 


The cantilever spring is placed so that it is unextended 
(zero force) when the tip is just in contact with the un¬ 
deformed bubble. 

The derivative of the extended total entropy with re¬ 
spect to the tip position, at constant tip contact radius 
ri, and evaluated at the optimum profile z(r;ri,zo) is 
now required. 

Above, in deriving the Eular-Lagrange equation for 
the profile, the variation at the boundaries vanished, 
Sz{ri) = 5z{r^) = 0. In the present case, the variation 
at contact on the tip must be dz{ri) = Azq. This means 
that one picks up an extra term from the integration by 
parts of the variation in area. 


rrs 

5A\z] = 27r / drr 
Jr^ 


z'ir) 


= —2'Kr 


ri \/l + z'(r) 

z'(ri) 


-Jz'(r) 


V^l -f z'(ri)2 


Azo 


( 2 . 20 ) 
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Accordingly 

^Siot,k 


T 


dzn 


’(t) 


A, 


SV[z] 6A[z] 

' Sz(r) 6z(r) 


+ ITT Ap + 


?(r-) 
27rriz'(ri) 
\/l + z'(ri )2 


dz{r) 

dzn 


- kt[zo - Zc] 

= 0 + F + Ft. 


( 2 . 21 ) 


The first term vanishes for the optimum profile. The re¬ 
maining term proportional to Ap arises from the deriva¬ 
tive of the constant contributions to the volume. This 
and the remaining term proportional to 7 give the force 
exerted by the bubble on the tip, 


F = TTT^Ap -I- 


27rri2:'(ri) 

V'l -b z'(ri)^^ 


( 2 . 22 ) 


This is just the pressure difference times the cross-section 
contact area plus the vertical component of the surface 
tension force times the contact perimeter. 

The quantity T) = —kt[zT) — Zc] is the force exerted by 
the tip on the bubble. One sees that in the equilibrium 
or static case, when the extended total entropy is a max¬ 
imum, the force due to the bubble is equal and opposite 
to the force due to the tip, F = — Ft(zo)- A positive bub¬ 
ble force, as sketched in Fig.[Tl corresponds to a negative 
cantilever force (in a signed sense). In practice, one often 
calls the cantilever force the applied force, or the load. A 
positive cantilever force (negative bubble force) gives an 
extended rather than a flattened bubble. 

Comparing bubble force with the first integral, 
Eq. (12.1811 . one sees that the integration constant is just 
the force exerted by the bubble on the tip, so that one 
has 


F = Trr^Ap -|- 


2TTrz'{r) 

^l + z'{rY^ 


(2.23) 


Hence one has a differential equation for the profile as a 
function of the force, F, the substrate contact radius ra, 
the tip contact radius ri, and the tip contact height zi. 
Given the geometry of the tip (eg. the half angle a), the 
location of the tip of the tip, zq, is determined by the 
latter two quantities, if it is ever required. 

In what follows an analytic expression for the profile 
will be derived in the weak force limit. The expansion is 
valid when \F\Rc/2TT^{s)r‘l <C 1 . 


2. Location of Dimple Rim 


The dimple rim r 2 is the maximum height of the bub¬ 
ble, so that z'{r 2 ) = 0. From the differential equation for 
the profile this yields 



(2.24) 


Note that since r 2 > ri, this sets a lower limit on the 
repulsive force that gives rise to a dimple. Of course there 
is also no dimple for attractive forces (negative loads). In 
either case, the dimple rim plays no further role in the 
analysis. 


3. Profile 


Rearranging equation (12.2311 for the profile, gives 
[F — 7rApr^]^[l -t- z'^r)"^] = (2TT'jr)‘^z'(r)'^, (2.25) 


or 


±[F - TrApr^] 
(27r7r)^ — [F — TrApr^]^ 


(2.26) 


The positive root is the physical root. 

One can see that a minimum and a maximum force is 
defined for a given pinned tip contact line ri when the 
gradient of the profile becomes infinite, z'(ri) = ± 00 . 
The bubble will rupture when the applied load exceeds 
these limits. From the profile equation one has 


and 


Fmin = -27rri7 


Fmax = 27rri7 


I - 


IL 

Rc 




(2.27) 


(2.28) 


In fact, since the gradient of the profile can’t be infinite 
anywhere, these two limits hold for any r on the interval 
[ri, ra]. Using ri gives the tightest upper bound because 
r > ri. But it can be the case that a tighter lower bound 
can occur by taking r inside the interval. In particular, 
if ri < (Fc/2) < Ta, then the bubble will rupture if F < 
— 7 ri?c 7 / 2 . If ra < Fc/2, then the bubble will rupture if 
F < -27rra7[l - ra/Fc]- 

For small loads, |F| <C 2'iT^(s)r\/ one can expand 
the profile equation to linear order in the force. 


z'{r) = 


F — TrApr^ 


Y^(27r7r)^ — (TrApr^)^ -p 2(7rApr2)F -|- 0{F'^) 


iF/nAp)-r^ 


r^Rl - r2 + 2(F/7rAp) 


—r 


F/nAp 


+ 


rF/irAp 


^Rff.—r'^ r^Rf^ — [Fc ~ 

F0{F'^). (2.29) 


Recall that Ap = 27 /Fc. The first term gives the un¬ 
deformed profile, Eq. (1^ . and the remainder give the 
perturbation due to the force to linear order. With 
e(r) = z{r) — Zc{r), this gives the derivative of the per¬ 
turbation to linear order, 

f( ] - FRc/‘2Trj rFFc/27r7 
~ ryF 2 ^+ [F2-r2]3/2' 


r2 = 


(2.30) 
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FIG. 2: Pinned bubble profiles for various applied loads. The 
dotted curve is the undeformed bubble (zero load), a full curve 
is an exact profile, and an adjacent dashed curve is the cor¬ 
responding linear approximation. Above and below the un¬ 
deformed profile the loads are ±0.5 nN, ±lnN, and ±2nN. 
The parameters are s = 4, j(s) = 0.018 N/m (s* = 5), and 
ra = 100 nm, giving Rc = 120 nm, and Zc = 54 nm. Each 
curve terminates at its last stable contact radius. 


The integral of this is 


e(r) 


_LJ1i i-\/r^ 1 1 

27r7'I 2 1 ± \/l — \/l — ’ 

(2.31) 


where x = r/Rc- The integration constant is determined 
by the condition that e(r 3 ) = 0, 


llnl- 4 . 

2 1 ± v^l-x§ 


(2.32) 


FIG. 3: Ratio of exact effective bubble spring constant, 
AFj/S.z{r\), to linear bubble spring constant fcb, Eq. (12.371) . 
as a function of the applied load. The solid, dashed, and dot¬ 
ted curves are for ri = 10, 20, and 30 nm, respectively. Other 
parameters as in Fig. (2] Each curve terminates at its limits 
of solution. The inset magnifies the region around zero force. 


laboratory frame, the position of the tip, the tip contact 
circle, and the undeformed bubble interface are 

Co ~ -^0 4” Cs; Cl ~ 4” Cs? and Cc — ± Csi (2.33) 

respectively. Initially, the substrate is at Cs = 0, the 
tip is just touching the undeformed bubble Co = Zc, and 
the spring attached to the tip is undeflected. Hence in 
general the deflection is i5t = Co ~ and force exerted 
by the tip on the bubble is 

Ft =-hSt =-h[Co - z,], (2.34) 

where kt is the tip spring constant. 

Now zi = Zc{ri) + s{ri), or 


Figure H] shows several profiles of deformed bubbles. 
The exact profile was obtained by numerical integra¬ 
tion of the profile equation, Eq. (12.261) . Here and below 
the linear approximation refers to analytic results based 
the expansion to linear order in the force, in this case 
Eq. (12.311) . In the case of the figure, 2t: r'l'y/= 0.4nN 
for ri = 20 nm. For loads with magnitude much less than 
this the linear approximation can be guaranteed accurate 
accurate. It can be seen in the figure that the perfor¬ 
mance of the linear approximation is rather better than 
is indicated by this parameter. For larger loads there is a 
significant discrepancy between the exact and the linear 
profile. The problem is more acute for extensive than for 
compressive forces. 


4- Bubble Spring Constant 

The various quantities 2 above were measured relative 
to the substrate. Now, in the laboratory frame of refer¬ 
ence, let Cs be the position of the solid substrate. In this 


Cl = CcCr-i) ±e(ri). (2.35) 

The amount of bubble deformation at contact is 

<ri) = (2-36) 

where F is the force exerted by the bubble. The bubble 
spring constant is given by the profile equation evaluated 
at ri, 


fcbub = -2ttj - In 


1 ± \/i — xf \/T 

This depends upon the two pinning i 

Rc > rs > ri. 



(2.37) 

Obviously 


Figure [3] shows the exact effective spring constant of 
the bubble, AF/Az{ri) obtained from the numerical in¬ 
tegration of the profile equation, Eq. (12.261) . normalized 
by the analytic expression obtained from the expansion 
to linear order in the force, Eq. ()2.37l) . The expansion 
is valid when |F| <C 2TT^{s)r\/Rc. For the conditions in 
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FIG. 4: Bubble linear spring constant as a function of tip 
contact radius xi = rijRc for a substrate radius *3 = rs/Rc 
of 0.9 (solid curve), 0.7 (dashed curve), and 0.5 (dotted curve). 
In each case, the upper limit is a;™®’' = [1 — -^1 — 2 : 3 ] tana, 
with a = 10 °. 


Fig. 131 the right hand side is 0.09, 0.38, and 0.85 nN for 
ri = 10, 20, and 30 nm, respectively. One can indeed see 
in the figure that the bubble behaves linearly to within 
about 20 % of the exact value when the loads lie within 
this bound. 

Figure [3] shows the ratio of the bubble linear spring 
constant to the surface tension. It can be seen that the 
variation is rather weak over the practical range, with the 
bubble spring constant being not more than a factor of 
two larger than the surface tension. In general the bubble 
spring constant is larger than the surface tension unless 
the contact radius is very small. Clearly, in the absence 
of specific information about the size of the pinning radii, 
one will not go too far wrong in taking the liquid-vapor 
surface tension of the supersaturated interface to be 0.5- 
1 times the measured bubble spring constant. 

In the linear regime, the bubble and tip act as two 
springs in series. To see this explicitly, consider a change 
in the position of the substrate, A^g, at constant ri (and 
ra). Since Ari = 0, the change in contact position must 
equal the change in tip position, A^i = ACq. In the static 
situation AF) = —AF, and so one has 


fctACo = -fcbubA£(ri) 

= -fcbub [ACi - ACc(ri)] 

= — fcbub [ACo — ACs] . 

This implies that [kt -I- fcbub] ACo = fcbubACs, or 

ACo _ fcbub 
ACs “t” ^bub 


(2.38) 


(2.39) 


This ratio is less than unity, whereas for a hard surface it 
would be unity. A measurement of the slope of the deflec¬ 
tion versus the drive distance allows the spring constant 
of the interface to be determined. If the contact radii 
are known, then the surface tension of the bubble can be 
estimated. 



FIG. 5: Hemispherical bubble penetrated and deformed by a 
conical hydrophobic tip (blunt radius ro) upon which it slips. 


Alternatively, in terms of the separation between the 
tip and substrate, zq = Co “ Csj one has 


ACo ^ ACo 

Azo ACo — ACs 

1 

“ 1 - ACs/ACo 

1 

^ _ fct~l~fcbub 

^bub 

_ ^bub 

kt 


(2.40) 


In words, the slope of the deflection versus separation 
curve equals the negative of the ratio of the bubble and 
cantilever spring constants. 


III. DEFORMED BUBBLE: PRICK SLIP 


For the case of slip on the tip, the liquid-vapor-solid 
contact circle is not pinned at ri (see Fig. [5|). One now 
has also to maximize the total entropy with respect to 
ri, taking into account the difference in solid surface en¬ 
ergies, Ay = 7 sg — 7 si. This is negative for a hydrophobic 
tip. 

For a macroscopic bubble or droplet on a planar sub¬ 
strate made of the same material as the tip, the equilib¬ 
rium condition is that the contact angle measured in the 
liquid phase satisfies 7 cos 6t = Ay. This does not hold in 
the present case of a conical tip penetrating the bubble. 

For the present case of slip, the total entropy is 


Ftot[ 2 ] = 5b(fV, U,T) - 1a[z] - |u[z] + 




(3.1) 


The surface area of the blunt conical tip inside the bubble 
is 




TTr\ 

sin a 


sin a 


-bTrrg. 


(3.2) 
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The formulae for the volume V[z] and area A[z] were 
given above and depend upon ri. Also recall that ri — 
fo = [zi — zq] tana. 


Differentiating with respect to the profile (at constant 
Ti) gives the Eular-Lagrange equation for the optimum 
profile, as given above. Differentiating with respect to 
ri, one has dri = dzitana, which means that Sz(ri) = 
dri/tana, and so the extra boundary term in the profile 
derivative of the area appears, as in 8 IIB II Also, dzo = 
0. In view of this one has 


dSu 


dri 


z{r) 

Ap 5V[z\ 7 5A\z] 


(3.3) 


dz{r) 

_ T 5z{r) T 5z{r)\-^^^ dri 


+ 


2'Kr 

tana 


- + 27rri (zo - 
a \ tana/ 


T 


-y /1 + z'{ri)^ tana 


27rr] 

Ay 

T sin a 

27rri j 

2 zi 

r 1 

Rc 


n - 77 — 7T zVri)/tana 


A 7 


sma 


7 

T 


Explicitly the force exerted by the bubble on the tip is 
,d5tot 


F = T- 


dzo 


= T- 




Sz{r) 




dz{r) 

dzn 


dStot{[z],ri,zo) 


dri 

dStot{[z],ri,zo) 


dri 


= T 


dzo 

dStoti[z],ri,zo) 


dzo 


= Trr^Ap + 


Zj,Zi 

27rrijz'(ri) 

+ z'(riy 


(3.6) 


Hence the expression for the profile is unchanged from the 
stick case, although of course for a non-stick hydrophobic 
tip, E < 0 . 


5. Algorithms 

Recall the equations of for the undeformed profile, 
Eq. (|2.5|) et seq. Also recall that in the linear approxima¬ 
tion the deformation e(r) = z(r) — Zc{r), and its deriva¬ 
tive s'{r) = z'(r) — z'(r), are linearly proportional to the 
force, Eqs (12.301) and (12.311) . In view of these one can 
succinctly write 


z(ri) = Zc(ri) - fcbub(D) (3.7) 


Setting this to zero, the trivial solution is ri = 0. The 
non-trivial solution gives an equation that the profile 
slope at the optimum radius must satisfy. 


Ay 

7 sin a 


2zi 


\/l + z'(ri)2 -I- 


z'(ri)/ tana 
^l + z'iri)"^ 


(3.4) 


This replaces the planar contact angle condition. 
In the limit a —^ 0 this is 
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and 

z'iri) = z '( e ) + qhuh{ri)F. ( 3 . 8 ) 

Write ki = khuhin) and qi = gbub(ri). 

To linear order the optimum contact radius satisfies 
(recall that ri — ro = [z{ri) — zo] tana) 
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which is the expected contact angle condition for a cylin¬ 
drical tip. 

The expression for the force ought to be unchanged. 


Rearranging this gives explicitly E(ri). 

Let the equilibrium curve be E(ri), Cs(^i)- The deflec¬ 
tion is (5t = Zo + Cs — Zc = — Et/fct = E/fct. The equilib¬ 
rium curve in linear approximation can be generated as 
follows: 









































































FIG. 6: Equilibrium (slip) deflection versus separation curves 
for a tip penetrating a nanobubble for different tip radii. 
These are the exact theory for, from bottom to top at contact, 
blunt tip radii of ro = 0, 20, 30, 40, and 50 nm. The dotted 
line is a guide to the eye. The tip has a = 10°, fct = 0.35 N/m, 
and Ay = ON/m or = 90°. All other parameters as in 
Fig.H 

• choose ri 

• calculate F{ri), and zq = z{ri) + [rp — rij/tana 

• calculate from F{ri) = —Ft = kt[zo + Cs — Zc] 

• plot deflection, (5t = F{ri)/kt versus separation, zq. 

For the exact, non-linear calculations, one has to spec¬ 
ify the load F and calculate the profile z{r-, F) by numer¬ 
ical integration from (r3,0). For the equilibrium (slip) 
case, one terminates the profile at the value of ri where 
the profile has the equilibrium contact angle, Eq. (laa . 
In the cases where there are two solutions one chooses 
the one based on continuity. From the value of ^(r’l), 
one obtains zq and Q as in the linear case. One then 
chooses a new load and repeats the process. For the case 
of the pinned tip contact line, for each F one instead 
terminates the profile at the fixed value of ri. 


6. Results 

Figure | 6 ] shows several equilibrium deflection versus 
separation curves. Equilibrium here and below mean that 
the contact circle on the tip is free to move. Hence the 
optimum contact angle that maximizes the entropy is 
established at each separation. Eor brevity this is also 
called slip. The data in the figure is obtained with the 
exact theory. Results obtained with the linear theory are 
entirely obscured by the exact curves on the scale of the 
figure. 

In the case of Eig. [ 6 l the tip has been taken to be 
indifferent to water. Ay = ON/m or = 90°, and for 
the most part the force is repulsive. This corresponds 
to a positive cantilever deflection and to a compressed 



FIG. 7: Equilibrium (slip) deflection versus separation for 
a blunt tip penetrating a nanobubble for different tip surface 
energies. The solid curves are the exact theory and the dashed 
curves are the linear theory (obscured in the top three cases), 
with ro = 20 nm and, from bottom to top. Ay = —10.0 mN/m 
{9t = 98°), -6.3mN/m {dt = 95°), -2.5mN/m [Ot = 92°), 
OmN/m {6t = 90°), and +2.5 mN/m {9t = 88 °). No exact 
solution was found for the lowest curve. The dash-dotted 
curve and the partially obscured dash-double dotted curve 
are the exact and linear results respectively for ro = 50 nm, 
and Ay = —6.3mN/m {6t = 95°). All other parameters as in 
Fig.H 


bubble, as in Fig. [T] and in the lower half of Fig. [5] The 
exception is for for the infinitely sharp tip, rg = 0 nm, 
which shows a weak attraction. 

It is noticeable that the force curves are almost straight 
lines, and that their magnitude increases with increasing 
tip radius. The predominant reason that the force in¬ 
creases with decreasing separation is the cone half angle, 
which means that the contact radius increases as the tip 
penetrates further into the bubble with decreasing sep¬ 
aration. The main reason that the force becomes more 
repulsive with increasing end radius rg is that the repul¬ 
sive pressure contribution is proportional to whereas 
the often attractive surface tension contribution contains 
a factor ri. 

Figure [7] explores the effect of the tip surface energy 
difference on the force curves. The conversion of the sur¬ 
face energy difference to a macroscopic tip contact an¬ 
gle uses the Young equation, the saturation value of the 
surface tension, y'l' = 0.072 N/m, and assumes that the 
difference in surface energies is unchanged by the level of 
supersaturation of the solution. Ay = yl cos Ot. 

It can be seen that as the surface energy difference 
becomes more negative, the force at a given separation 
becomes increasingly attractive. The penetrated bubble 
is extended from its undeformed shape (see Fig. [S] and 
the upper half of Fig. HD. The slope of the almost linear 
deflection versus separation curves changes from nega¬ 
tive to positive for Ay < —5.0 mN/m {9t > 92°), close to 
where the surface energy difference changes sign, and it 
becomes increasingly positive as the surface energy dif- 
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ference becomes increasingly negative. 

The linear theory becomes increasingly less accurate as 
the surface energy difference becomes more negative. In¬ 
deed, whilst the linear theory produced a solution curve 
for Ay < —lO.OmN/m at this tq = 20 nm, no exact solu¬ 
tion was found (essentially because the bubble ruptured 
at this contact radius). 

Figure [7] also presents the case of tq = 50 nm and 
Ay = —6.3mN/m {dt = 95°). Compared to the same 
surface energy but for rp = 20 nm, one sees that the 
force is less attractive and as well the slope has decreased, 
changing from positive to negative. As mentioned above, 
increasing the contact radius increases the repulsive pres¬ 
sure contribution to the force more than the attractive 
surface tension contribution. One can conclude that a 
negative slope is the signature of pressure dominance 
(large contact radius, small magnitude or positive surface 
energy difference) whereas a positive slope is the signa¬ 
ture of surface energy dominance, (small contact radius, 
large negative surface energy difference). 

In Figs inland [7] it can be seen that there is a discontin¬ 
uous jump from zero deflection prior to contact with the 
undeformed bubble to a non-zero deflection immediately 
after contact. This is due to the fact that non-contact 
forces are neglected in the present model and also to the 
fact that the end of the tip has been taken to be a disc 
(ie. perfectly blunt, planar). In the present calculations 
liquid-vapor interface contact with the flattened end has 
been excluded. The size of this jump increases with the 
tip end radius rp. Obviously this is an idealized model 
of the actual tip of the tip, which is actually neither per¬ 
fectly sharp nor perfectly blunt, and in practice there 
may be a smooth transition rather than a jump. The 
cantilever manufacturer typically quotes a tip radius 20- 
60 nm. Also of course the AFM cantilever tip is generally 
in the form of a square pyramid rather than the present 
right circular cone. 

Figure [8] shows both a slip trajectory and a stick-slip 
trajectory for a blunt tip, rp = 20 nm. The stick branches 
were chosen more or less randomly, with one eye on aes¬ 
thetics, one eye on fundamental considerations, and one 
eye on experimental data. There is no fundamental rea¬ 
son that when the contact line gives way it should jump 
to, and immediately stick at, the equilibrium position. It 
could stick prior to reaching the equilibrium position, or 
it could slip along the equilibrium curve after the jump. 
One could perhaps argue for a yield stress such that the 
contact line always slipped when the excess force per unit 
contact line reached a certain value, but this has not been 
done here. 

The first tip contact line stick has been taken to occur 
at ri = 26.4 nm. For the exact calculation, the initial 
slope of the stuck branch is -0.078, and the average slope 
of the entire branch is -0.090. The linear prediction for 
the slope at this contact radius and surface tension is 
-0.137. 

The second tip contact line stick occurs at ri = 
28.9 nm, and the exact initial slope is -0.86, and the av- 
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FIG. 8: Deflection versus separation for a tip penetrating 
a nanobubble. The solid curve is the equilibrium slip case 
using the exact theory, and the parallel dashed curve uses the 
linear theory. The dash-dotted curve is a stick-slip case using 
the exact theory, with the pinned contact radii being ri = 
26.4, 28.9, and 31.6 nm (chosen more or less arbitrarily). The 
approach curve overlaps the initial part of the retraction curve 
on the final stick region. The dotted arrows show cantilever 
jumps. The final jump out is at the limit of bridging bubble 
stability for ri = 31.6 nm. The blunt tip radius is ro = 20 nm, 
and the difference in tip surface energies is Ay = —6.3mN/m 
{9t = 95°), with all other parameters Fig. [S] 


erage slope of the entire branch is -0.95. The linear pre¬ 
diction for the slope at this contact radius is -0.143. 

The third tip contact line stick occurs at ri = 31.6 nm, 
and the exact initial slope of the stuck branch is -0.083, 
and the average slope (extension, into contact) is -0.093. 
The linear prediction for the slope at this contact radius 
is -0.150. 

For the single stuck branch on retraction, ri = 31.6 nm, 
curvature is quite evident, which graphically indicates 
that the linear theory is inapplicable. This curvature, 
including the final flat region, is consistent with the data 
for the effective exact bubble spring constant in Fig. [S] 
The vanishing slope is due to the fact that the bubble 
spring constant tends to zero with increasing bubble ex¬ 
tension. 

It is evident that in the case treated in Fig. [U the ra¬ 
tio of the exact slope to the linear slope is about 2:3. 
This is in agreement with the ratio indicated in Fig. [3] 
for F = —1 nN. Fitting the linear approximation to mea¬ 
sured experimental data where stick is evident, partic¬ 
ularly on the extension branch close to the equilibrium 
curve, can be expected to provide a useful first estimate 
for the bubble surface tension. 
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IV. EXPERIMENT 
A. Limitations of Theory 

Before analyzing any experimental data, it is worth 
enumerating the limitations of the present theory. 

First, the theory models the tip as a right circular cone, 
whereas in reality an AFM tip is a rectangular pyramid. 
Second, it models the tip of the tip as perfectly blunt, a 
disc of radius ro, whereas in reality the tip may be curved 
from wear. Third, it is assumed that the cross-section of 
the nanobubble is perfectly circular, which need not be 
the case if substrate heterogeneity determines where the 
contact line is pinned. Fourth, it is assumed that the sub¬ 
strate contact line does not alter during the penetration 
of the nanobubble by the tip. 

Fifth, it is assumed that the tip penetrates the 
nanobubble at the apex, whereas in reality it can pen¬ 
etrate off-axis, either by design or by accident. Sixth, it 
is assumed that the tip is oriented normal to the sub¬ 
strate, whereas in reality the cantilever and tip are tilted 
at about 11° in the axial plane. Seventh, it is assumed 
that the air inside the nanobubble remains in diffusive 
equilibrium with that in the solution during the force 
measurement (constant chemical potential), whereas in 
reality the measurement might be rapid enough for con¬ 
stant number to hold instead. 

The fifth and sixth points mean that the predicted nor¬ 
mal force on the cantilever has an error that increases 
with displacement from the central axis of the nanobub¬ 
ble, that the displacement from the central axis varies 
with the separation and with the deflection of the can¬ 
tilever during a force measurement (due to the tilt), 
and that there is a torque on the cantilever due to the 
asymmetric forces on the tip when it passes through the 
nanobubble interface off-axis. This nanobubble torque is 
in addition to the torque that acts in all AFM force mea¬ 
surements on the tip in contact with the hard substrate 
due to the normal and lateral (friction) forces, which are 
implicitly included in the photo-diode calibration.— “— 

The limitation summarized in the fifth point can be 
substantially alleviated by ensuring that the separation 
at which first contact with the nanobubble occurs is equal 
to the measured height of the undeformed nanobubble. 
One is aided in this by the fact that the nanobubble pro¬ 
file is horizontal near the apex. In this case one can be 
confident that, at least for small cantilever deflections 
and separations close to initial contact, the tip is pene¬ 
trating the nanobubble close to the apex, and the cali¬ 
bration factor is correct. In this regime the present cal¬ 
culation of the nanobubble normal force and the neglect 
of any nanobubble lateral forces or torques, ought to be 
accurate. The consequences of the seventh point are dis¬ 
cussed below on p. El 

In view of these limitations, one ought not to expect 
the present theory to be able to quantitatively describe 
every aspect of an individual nanobubble force measure¬ 
ment. The main goal in the first instance is to establish 


that the nanobubble surface tension is less than the usual 
air-water surface tension, and if possible to quantify re¬ 
liably its value for a given solution. 

To this end the following protocol was adopted. Pri¬ 
mary emphasis was placed on the first pinned region of 
the measured separation-deflection curve, since this is the 
one that can be guaranteed closest to the apex. Addi¬ 
tional pinned regions were used to confirm the values 
deduced from the first one. 

Since the surface tension that is required to fit a given 
slope decreases with increasing value of ri ^ ro, one can 
establish an upper bound for the surface tension by spec¬ 
ifying the lowest realistic value of tq . In view of the spec¬ 
ification that a new tip has radius in the range 20-60 nm, 
fixing for example tq = 20 nm should give an upper limit 
on the surface tension. 

Further, as is discussed below on p. 1111 for a given 
surface tension and pinned contact radius ri, the slope 
calculated at constant chemical potential is less in magni¬ 
tude than the slope calculated at constant number. This 
means that the surface tensions obtained below at con¬ 
stant chemical potential are larger than those that would 
be required to ht the slopes at constant number. Again 
one can be confident that the surface tensions obtained 
here are an upper bound on those for actual nanobubbles. 

Since the slope of the pinned regions equals the neg¬ 
ative of the ratio of the nanobubble spring constant to 
the cantilever spring constant, the surface tension can 
now be determined using the linear theory. The accu¬ 
racy of this can be checked against the exact theory. 
From the surface tension and the nanobubble curvature 
radius determined by tapping mode imaging, the super¬ 
saturation ratio is now determined, Eq. m- Using 
a linear model for the supersaturated surface tension, 
7 (s) = (s^ — s)7’l'/(s^ ~ i) ^20i26i27 = 0.072 N/m 

is the saturated surface tension, the spinodal supersatu¬ 
ration ratio can now be determined. 


B. Nanobubble 1 

Figure [9] shows AFM measurements of the force on 
a cantilever tip due to a single nanobubble. What is 
plotted is the positional deflection of the cantilever; to 
obtain the force, multiply by the cantilever spring con¬ 
stant, kt = 0.35 N/m. The first measurably significant 
deflection occurs at a separation of zq = 32.9 nm. This 
is in good agreement with the height of this particular 
nanobubble imaged in tapping mode, Zc = 33 nm. This 
suggests that pre-contact forces (van der Waals, electric 
double layer) are negligible. It also confirms that the 
measurement was performed in the central region of the 
nanobubble close to the apex. 

The nanobubble profile obtained from the image (not 
shown) has a contact radius of = 108 nm, which, with 
its height, corresponds to a contact angle of 146°. The 
fact that this contact angle is substantially higher than 
the contact angle of a macroscopic water drop on HOPG, 
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FIG. 9: Cantilever deflection versus separation for a SiN 
tip penetrating a nanobubble on an HOPG substrate on ap¬ 
proach (triangles) and retraction (crosses). The cantilever 
has spring constant kt = 0.35 N/m and conical half-angle 
a = 10°. The undeformed nanobubble has measured height 
Zc = 33 nm and substrate contact radius rs = 108 nm, corre¬ 
sponding to a curvature radius = 192 nm and a contact 
angle of 146°. The dotted lines and arrows are guides to the 
eye, with the adjacent number giving the slope. The full and 
dashed curves are respectively the calculated exact and linear 
equilibrium deflection (tip contact line slip) with ro = 10 nm, 
7 = 0.040 N/m, s = 5.16, and. Ay = — lO.OmN/m (equiva¬ 
lently, Ot = 98°). The dash-double dotted line (ri = 10 nm, 
ro = 6 nm, obscured), and the dash-dotted curve (ri = 20 nm, 
ro = 10 nm), are calculated exact deflections with the tip con¬ 
tact line pinned, using 7 = 0.044 N/m, s = 5.55. 

64-92° )^^!^^ is strong evidence that the nanobubble con¬ 
tact rim is pinnedi^ 

In the experimental data just after first nanobubble 
contact, the positively sloped region, followed by a brief 
plateau, followed by a small jump to the base of the first 
marked linear region, are all due to the initial spreading 
of the nanobubble on the tip of the tip and up its sides. 
This region is not well-modeled by the present geometry 
of the tip of the tip as a perfectly planar circular disc. 
The fact that the deflection is negative in this region 
indicates that it is favorable for the tip to penetrate the 
nanobubble, which is to say that the SiNi tip must be 
hydrophobic, or possibly barely hydrophilic. 

The two linear regions with labeled slopes evident in 
the experimental deflection data on approach confirm 
that the bubble can have a pinned contact line and be¬ 
have as a Hookean spring. Since the cone half angle is 
small, a = 10 °, to leading order one can take the contact 
radius used in fitting the slopes to be the same as the 
radius of the perfectly blunt tip, ri ~ tq. The manufac¬ 
turer’s quoted radius of the tip, 20-60 nm, which might 
refer to either the tip’s width or else its radius of curva¬ 
ture, can be assumed to be of the same order as rg in the 
present simple model. 

Choosing a contact radius at the lower end of the re¬ 
alistic values, ri = 10 nm, a value of 7 = 0.040 N/m 
gives k^i/kt = 0.27, which is the negative of the measured 


slope of the first linear region. Using Eq. this sur¬ 

face tension requires a supersaturation ratio of s = 5.16 
to give a critical radius of Rc — 193 nm, equal to that 
deduced from the tapping mode images of this partic¬ 
ular nanobubble. Conversely, choosing the upper limit 
ri = 50 nm, the fitted surface tension is 7 = 0.015 N/m 
and s = 2.58. Hence even in the most pessimistic case 
of smallest contact radius one can see that the nanobub¬ 
ble surface tension, 7 = 0.040 N/m almost a factor of 
two smaller than the surface tension of the saturated air- 
water interface, 7 = 0.072 N/m, and that the solution is 
substantially supersaturated with air, s = 5. 

The just quoted slopes were obtained with the linear 
theory. Applying the exact non-linear theory with the 
contact line again pinned at ri = 10 nm, the tangent 
at the start of the dash-double dotted curve in Fig. [51 
gives the required slope —0.27 using 7 = 0.044 N/m and 
s = 5.55. In this case using rg = 6 nm shifts the curve 
laterally to coincide with the measured data. Obviously 
using larger values of ri will require smaller values of 7 ( 5 ) 
to fit the slope. There is very little curvature evident in 
the non-linear curve. The 10% difference in the surface 
tension between the exact and the linear fits means that 
the linear theory provides an acceptable estimate of the 
surface tension from the slope that is both analytic and 
reliable. 

The slope of the second linear region in Fig. [51 —0.29, 
is fitted by 7 = 0.043 N/m using ri = 10 nm using the 
linear theory. Alternatively, the change in contact posi¬ 
tion of the nanobubble on the tip may be approximated 
as the change in separation at the base of the two lin¬ 
ear regions, Azi « Azg = 11 nm, assuming that the 
bubble profile is essentially the same in the two cases, 
which it would be if the contact line were mobile prior 
to the start of the pinned regions. The change in con¬ 
tact radius is Ari = Azi tan a « 2 nm. Using this one 
finds that a single surface tension (to better than 0.03%) 
7 = 0.0402 N/m gives a slope of -0.27 using n = 10 nm, 
and a slope of -0.29 using ri = 12 nm. In the case of 
ri = 50 nm and ri = 52 nm, the two slopes are given by 
a single surface tension with a slightly worse variance of 
2%. This tends to suggest that the smaller contact radius 
is more applicable, but this is by no means conclusive. In 
any case, that a single surface tension combined with the 
geometry of the tip fits the two slopes supports the model 
and the value of the surface tension. 

Here and below the calculations are performed at con¬ 
stant chemical potential, which is computationally conve¬ 
nient. Although there is strong evidence (see below) that 
the nanobubble is in diffusive equilibrium with the solu¬ 
tion over the time of the series of force measurements, it 
is unclear whether it is best to model each force measure¬ 
ment as at constant chemical potential or as at constant 
number. For the purposes of comparison, some exact 
calculations have been carried out at constant number. 
Using a surface tension of 7 ( 5 ) = 0.044 N/m, at a pinned 
radius of ri = 10 nm the tangent at zero force at constant 
chemical potential is -0.30, compared to -0.41 at constant 
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number. Using instead ri = 20 nm and the same surface 
tension, the tangent at zero force at constant chemical 
potential is -0.40, compared to -0.62 at constant num¬ 
ber. One sees that the slope has a higher magnitude at 
constant number, and that it increases relatively more 
rapidly with contact radius. Hence one would require 
smaller surface tensions to fit the measured slopes if one 
used constant number. The calculations here and below 
are at constant chemical potential, and so the surface ten¬ 
sions obtained represent an upper bound on the actual 
nanobubble surface tension. 

The two almost horizontal curves in Fig. IH] are the cal¬ 
culated exact and linear equilibrium curves, which as¬ 
sume that the contact line is mobile on the tip. The good 
agreement between the linear and the exact calculations 
is somewhat better than that for the predicted effective 
bubble spring constant. Both equilibrium calculations 
use ri = 10 nm, 7 = 0.04 N/m, and s = 5.16. In ad¬ 
dition a surface energy difference of A 7 = — lO.OmN/m 
was fitted, which corresponds to a macroscopic contact 
of water on a planar SiNi substrate of dt = 98°. This 
is slightly hydrophobic. The criterion for the fit, which 
was done by eye, was that the curve should pass close to 
the base of the two linear regions. (For reasons that are 
discussed below, the end point of the final jump was not 
included in the fit.) Since the cantilever jumps to these 
bases, the nanobubble at contact must also be moving 
along the tip, and so it can probably be assumed that 
the contact position is the equilibrium one. Of course 
the contact line may become pinned at the end of the 
jump prior to achieving its equilibrium position, so this 
fit may underestimate the magnitude of the surface en¬ 
ergy difference. Also, using a larger contact radius would 
require a smaller in magnitude value for the fit. Fortu¬ 
nately, the surface tension obtained by fitting the slopes 
of the linear regions is not affected by the tip solid surface 
energies. 

From the fitted value of the surface tension at ri = 
10 nm, 7 = 0.040 N/m, and the measured nanobubble 
curvature radius, Rc = 193 nm, which is equal to the 
critical radius, Eq. the supersaturation ratio can 

be deduced to be s = 5.2. The linear model for the su¬ 
persaturated surface tension is 7 ( 5 ) = (s^ — (s^ — 1 ), 

where 7 ^ = 0.072 N/m is the saturated surface tension, 
and is the spinodal saturation ratio. This has been 
shown to fit the available computer simulation data rea¬ 
sonably accurately]These computer simulations 
give Ri4-6 for a Lennard-Jones fluid, depending on 
the temperature. The present fit, 7 ( 5 . 2 ) = 0.040 N/m, 
gives the spinodal supersaturation ratio = 10.4. Al¬ 
ternatively, at ri = 50 nm, the linearly fitted supersat¬ 
urated surface tension 7 = 0.015 N/m requires s = 2.6 
to give the measured nanobubble curvature radius, and 
corresponds to = 3.0 in the linear model. 

Also shown in Fig. |9] is the calculated exact (ie. non¬ 
linear) deflection on extension with the contact radius 
pinned at ri = 20 nm. The value of the contact ra¬ 
dius was chosen so that the calculated curve fitted by 


eye the measured data at the end of the retraction 
branch. This calculation used the non-linear fitted value 
7 = 0.044 N/m and s = 5.55, and also tq = 10.3 nm, 
which shifts the curve laterally. The consistency of this 
with the exact fit to the slope —0.27 (ri = 10 nm and 
ro = 6.2 nm) is probably already acceptable; with a little 
optimization of ri it could doubtless be made even bet¬ 
ter. It is also undoubtable that larger values of ri and rg 
and smaller values of 7 ( 5 ) could also fit both the slope of 
the pinned regions on approach and the flattened region 
on retraction. The conclusion that one can draw is that 
for separations zq > 20 nm the retraction data can be 
described by the pinned nanobubble model using param¬ 
eters consistent with what was deduced from the slopes 
of the pinned approach data. 

For separations zq ^ 20 nm on retraction, and zq ^ 
3nm on approach, the data in Fig. [9] is not described by 
either the pinned or the mobile contact line theory. Sim¬ 
ilar steep curved regions have been observed in a number 
of other AFM nanobubble measurements, albeit for col¬ 
loid probes rather for tipsi^i^ The origin of this particu¬ 
lar behavior is unclear. Because of the coincidence of ap¬ 
proach and retraction here, this is clearly an equilibrium, 
non-dissipative phenomenon. Calculations show that it 
is not due to elastic deformation of the substrate (not 
shown)It might be due to torque on the cantilever, 
although measurements across the nanobubble indicate 
that this is in general negligible. That the curve is much 
steeper at F = 0 than the clearly pinned regions suggests 
that it is not due to pinning of the contact line, unless 
the pinned contact radius had increased very substan¬ 
tially. The apparently contiguous flat region r ^ 20 nm 
on retraction is similar to the non-linear calculations of 
the force due to pinning of a highly extended nanobub¬ 
ble with small contact radius. Obviously whatever the 
origin of this behavior, it could be simply additive to the 
force due to the pinned (or slipping) contact line since 
the nanobubble force is always present. Because of the 
uncertainty as to the origin of this force close to contact 
it has been neglected in fitting the nanobubble. 

The measurements in Fig. [9] were part of a sequence 
of twelve successive force measurements across this par¬ 
ticular nanobubble (not shown). The number, position, 
and extent of the linear regions could differ between force 
measurements, presumably because contact line stick and 
slip are stochastic events, but the slopes were unchanged. 
This suggests that torque on the cantilever due to off- 
apex penetration has negligible effect. Tapping mode 
images before and after the sequence of force measure¬ 
ments show that the nanobubble itself was unchanged in 
size and shape by the force measurements. This is strong 
evidence that the nanobubble is thermodynamically sta¬ 
ble, and that even penetrating it a dozen times with the 
cantilever tip did not destroy or alter it. It is also evi¬ 
dence that the nanobubble is pinned at it contact rim. 
In several other series of measurements, up to a hundred 
force measurements were performed on a single nanobub¬ 
ble, interspersed with several AFM tapping mode images, 
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FIG. 10: Measured cantilever deflection versus separation for 
a nanobubble (curves and lines as in preceding figure; dif¬ 
ferent nanobubble, cantilever, and solution). The cantilever 
has spring constant fet = 0.24 N/m, and conical half-angle 
a = 10°. The undeformed nanobubble has measured height 
Zc = 9.6 nm and substrate contact radius ra = 82.5 nm, corre¬ 
sponding to a curvature radius = 359 nm and a contact an¬ 
gle of 167°. The calculated equilibrium exact (solid) and lin¬ 
ear (dashed, obscured) curve use ro = 10 nm, 7 = 0.041 N/m, 
s = 3.3, and A 7 = -1-1.3mN/m (equivalently, 9t = 89°). 

and no significant change in the nanobubble was observed 
in any case. 


C. Nanobubble 2 

Figure [TO] shows results for another nanobubble in a 
different solution, with stick-slip behavior evident. As¬ 
suming an initial tip contact radius of ri = 10 nm, the 
measured slope —0.50 corresponds to a surface tension of 
7 ( 5 ) = 0.041 N/m (linear approximation). Using instead 
ri = 20 nm gives 7 ( 5 ) = 0.028 N/m. Larger contact radii 
require even smaller surface tension to yield this slope. 
Exact calculations differ by less than 1% from the linear 
results for the bubble spring constant in this regime. 

The second pinned region with slope -0.57 corresponds 
to 7 ( 5 ) = 0.047 N/m using ri = 10 nm, and to 7 ( 5 ) = 
0.032 N/m using ri = 20 nm. The third pinned region 
with slope -0.66 corresponds to 7 ( 3 ) = 0.054 N/m using 
ri = 10 nm, and to 7 ( 5 ) = 0.037 N/m using ri = 20 nm. 

From the point at which the pinned regions extrapo¬ 
late to zero deflection, one can deduce the value of the 
pinned radius ri for a specified value of the tip radius 
tq. However, it is not possible to find a single value of 
To which yields a single surface tension when all three 
slopes are fitted. For example, fixing vq = 10 nm, one 
finds that the three slopes -0.50, -0.57, and -0.66 corre¬ 
spond to ri = 10.9, 11.1, and 11.7nm, and to 7 ( 5 ) = 
0.039, 0.045, and 0.050 N/m, respectively. There is noth¬ 
ing wrong with ri increasing with each successive pin¬ 
ning event as the tip penetrates the nanobubble, but one 
would have hoped for a single surface tension. The best 



FIG. 11: Measured cantilever deflection versus separation for 
a nanobubble (curves and lines as in preceding figure; dif¬ 
ferent cantilever, solution, and nanobubble). The cantilever 
has spring constant fet = 0.35 N/m, and conical half-angle 
a = 10°. The undeformed nanobubble has measured height 
Zc = 19.5 nm and substrate contact radius ra = 205 nm, cor¬ 
responding to a curvature radius Rc ~ 1087 nm and a contact 
angle of 169°. The dashed curve is equilibrium (slip) lin¬ 
ear theory with ro = 20 nm, 7 = 0.037 N/m, s = 1.68, and 
A 7 = -1-3.8 mN/m (equivalently, 9t = 87°). 

that can be done (ie. minimizing the sum of the rela¬ 
tive standard deviations in surface tension and in tip ra¬ 
dius), yields rg = 12.323 ± .001 nm and 7 ( 3 ) =0.036, 
0.041, and 0.046 N/m, respectively. Possibly doing the 
calculations at constant number rather than the present 
constant chemical potential might yield more consistent 
results (see p. mi). Or possibly the problem is related 
to the unknown origin of the steep hook at small sep¬ 
arations discussed on p. In any case, the preferred 
value of surface tension is the one taken from the very 
first pinned region, since this lies closest to the tip of the 
tip of the cantilever and to zero force. 

The surface tension obtained using ri = 10 nm, 7(3) = 
0.041 N/m, and the nanobubble radius of curvature Rc = 
359 nm correspond to a supersaturation value of 3 = 3.3. 
Using this in the linear model for the supersaturated 
surface tension gives a spinodal supersaturation ratio of 
3^ = 6.3. Using instead ri = 20 nm, 7(3) = 0.028 N/m, 
correspond to a supersaturation value of 3 = 2.6. and a 
spinodal supersaturation ratio of 3 ^ = 3.6. 

Figure [TO] also shows equilibrium calculations using 
ro = 10 nm, 7(3) = 0.041 N/m, and 3 = 3.3. The value of 
the surface energy difference, A7 = -1-1.3 mN/m (equiv¬ 
alently, 0t = 89°) was chosen so that the curves passed 
through the base of the first pinned region. The exact 
and the linear calculations are almost indistinguishable. 


D. Nanobubbles 3 and 4 

Figure [TT] shows yet another nanobubble measurement. 
The AFM fluid cell was flushed with ethanol and then 
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water which means that there was likely exothermic 
heating. The undeformed nanobubble has apex height 
Zc = 19.5 nm with the first indication of a force occur¬ 
ring at a separation of 21.9 nm. The nanobubble was re¬ 
imaged after 47 force measurements and the apex height 
was Zc = 21.5 nm, and the substrate contact radius was 
ra = 215 nm, which correspond to a curvature radius 
i?c = 1081 nm. Again this is unambiguous evidence that 
the nanobubble is thermodynamically stable. 

The two linear regions measured in the figure, one each 
on extension and retraction, may be attributed to stick. 
Assuming a contact radius of ri = 20 nm, the slope —0.28 
corresponds to 7 ( 5 ) = 0.037 N/m. Alternatively for this 
slope, ri = 10 nm corresponds to 7 ( 5 ) = 0.048 N/m. and 
ri = 50 nm 7 ( 5 ) = 0.022 N/m. The exact and linear 
bubble spring constants agree to better than 0.05% in 
this regime. 

The value 7 ( 5 ) = 0.037 N/m, and the nanobubble ra¬ 
dius of curvature Rc = 1070 nm correspond to a super¬ 
saturation value of s = 1.7. Using this in the linear 
model for the supersaturated surface tension gives a spin- 
odal supersaturation ratio of = 2.4. Alternatively, 
7 (s) = 0.048 N/m corresponds to s = 1.9 and = 3.6. 

One can assume that the cantilever initially jumps in to 
the equilibrium slip position, which is supported by the 
flat nature of the deflection curve and the coincidence 
of approach and retraction. Assuming again a flattened 
conical tip with rp = 20 nm, these data are well-fitted 
using 7 ( 5 ) = 0.037N/m, s = 1.7, and A 7 = -1-3.8mN/m, 
which would correspond to a macroscopic contact an¬ 
gle of water on silicon nitride oi 9 1 = 87°. The va¬ 
lidity of the fit is supported by the coincidence of the 
jump-out separation and the end of the stability of the 
extended nanobubble. It can be mentioned that a vir¬ 
tually identical equilibrium slip curve can be obtained 
for rp = 50nm, with 7 ( 5 ) = 0.022N/m, s = 1.41, and 
with A 7 = -1-2.3 mN/m, which would correspond to a 
macroscopic contact angle of water on silicon nitride of 
dt = 88.2°. The equilibrium data can also be fitted by 
rp = 10 nm, with 7 ( 5 ) = 0.048 N/m, s = 1.9, and with 
A 7 = -k2.5mN/m, (0* = 88 °). 

Figure [12] shows another nanobubble in the same solu¬ 
tion as Fig.[TT] Despite the differences in height and con¬ 
tact radii, the two nanobubble have about the same ra¬ 
dius of curvature {Rc = 1087 nm there and 970 nm here). 
This is consistent with the level of supersaturation of the 
solution being unchanged and with the nanobubbles be¬ 
ing in diffusive equilibrium. Likewise, one would expect 
the surface tension to be unchanged, and this is confirmed 
by the fact that the slope of the linear pinned regions are 
about the same (—0.28 there and —.29 here). 

Using a contact radius of rp = 20 nm, the slope of the 
stick region in extension in Fig. [T^] of —0.29 corresponds 
to a surface tension of 7 = 0.035 N/m, a supersaturation 
ratio of s = 1.72, and a spinodal supersaturation ratio of 
= 2.4. The minimum of the extension curve touches 
a calculated equilibrium slip curve for A 7 = -1-2.5 mN/m 
(equivalently, 9t = 88 °). 



FIG. 12: Measured cantilever deflection versus separation for 
a nanobubble (same cantilever and solution as in the preced¬ 
ing figure; different nanobubble). The undeformed nanobub¬ 
ble has measured height Zc = 15.9 nm and substrate con¬ 
tact radius ra = 174 nm, corresponding to a curvature ra¬ 
dius Rc = 970 nm and a contact angle of 170°. The dashed 
curve is equilibrium (slip) linear theory with rp = 20 nm, 
7 ( 5 ) = 0.035N/m, and s = 1.72, and A 7 = -|-2.5mN/m 
(equivalently, dt = 88°). 


Using instead a contact radius of rp = 50 nm, the slope 
of the stick region in extension in Fig. [12] of —0.29 cor¬ 
responds to 7 ( 5 ) = 0.020 N/m, s = 1.42, = 1.58. Fit¬ 

ting the minimum of the extension curve gives A 7 = 
-1-1.3mN/m (equivalently, 0t = 89°). 

Using instead a contact radius of rp = 10 nm, the slope 
of —0.29 corresponds to 7 ( 5 ) = 0.046 N/m, s = 1.95, 
= 3.61. Fitting the minimum of the extension curve 
gives A 7 = -1-3.8 mN/m (equivalently, 0t = 87°). 

It should be mentioned that the slope of the stick re¬ 
gion was checked against that given by the exact theory 
and the agreement was better than 1%. The origin of the 
non-linearity and peak in the putative pinned region in 
Fig. [T^ (and in Fig. [TT|) is unclear, although one could 
speculate that the contact line might be moving with a 
finite velocity in these regions. 

The slight negative slope in the putative equilibrium 
curve here in Fig. [12] is difficult to reproduce in the the¬ 
oretical calculations. 


V. CONCLUSION 

The present paper gives analytic expressions for the 
nanobubble spring constant that allow its surface tension 
to be obtained from the slope of the pinned regions in 
a force-separation AFM measurement. Expressions are 
also given that allow the difference in tip surface energies 
(tip contact angle) to be obtained from an equilibrium 
part of the force curve. 

The present fits to the experimental data do not give 
enough information to pin down the value of the blunt 
tip radius rg. However, sensible results are generated 
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TABLE I: Measured and Deduced Properties of Nanobubbles 
and Solutions (ro = 10 nm). 


Figure 

kt 

(N/m 

rs Zc 

) (nm) (nm) 

Rc 

(nm) 

7(s) 

(N/m) 

S 


6»t 

(deg.) 

m 

0.35 

108 

33 

192 

0.040 

5.2 

10.4 

98 

[IQ] 

0.24 

82.5 

9.6 

359 

0.041 

3.3 

6.3 

89 

Er 

0.35 

205 

19.5 

1087 

0.047 

1.9 

3.6 

87 

M 

0.35 

174 

15.9 

970 

0.046 

2.0 

3.6 

88 


* Same solution, different nanobubble. 


by assuming a tip and contact radius ro = ri = 10 nm 
(Table [T|. This is at the lower end of the range of a typical 
AFM tapping mode cantilever tip, 20-60 nm. A high 
value of the radius, ro = ri = 50 nm gave quite low 
values of the surface tension, supersaturation ratio, and 
spinodal supersaturation ratio. 

The most reliable data appears to be that of Fig. [HI 
In this case the two pinned regions both appear to begin 
from the equilibrium curve, which means that the change 
in contact radius can be found from the change in separa¬ 
tion of the starts. Hence one has three knowns (the two 
slopes and the change in contact radius) and three un¬ 
knowns (the surface tension and the two contact radii). 
Solving this system yields the contact radii ri = 9.4 nm 
and r'l = 11.4 nm, and also 7 ( 3 ) = 0.041 N/m, s = 5.2 
and = 10.8. Unfortunately the other figures do not 
have pinned regions starting from the equilibrium curve 
and so the change in contact radius cannot be readily 
deduced. 

The estimates of the surface energy difference and 
macroscopic tip contact angle are based on the equilib¬ 
rium curves and are not so reliable. Little more can be 
said than that the contact angle is close to 90°. 

Likewise the estimate of the value of the spinodal su¬ 
persaturation ratio has limited reliability because of the 
simplicity of the linear supersaturated surface tension 
model. It is nonetheless consoling that it comes out to 
be of the same order as has been found in computer sim¬ 
ulations of a Lennard-Jones fluid.— 


One of the purposes of this study was to establish 
experimentally that the surface tension of nanobubbles 
was less than that of saturated water. The surface ten¬ 
sion is obtained from the slope of the pinned regions and 
does not rely upon the surface energy difference nor the 
spinodal supersaturation ratio. The greatest uncertainty 
concerns the tip radius, and to this end it is better to 
use a low value, since this overestimates the surface ten¬ 
sion. (Doing the calculations at constant chemical po¬ 
tential rather than at constant number further overesti¬ 
mates the surface tension.) Hence the data in Table U 
for ro = 10 nm, which give 0.04 < 7 ( 3 ) ^ 0.05 N/m, are 
most likely an upper bound on the nanobubble surface 
tension. The nanobubble surface tension is substantially 
reduced from that of the saturated air water interface, 
7 ^ = 0.072 N/m. 

The solution supersaturation ratios are deduced to be 
in the range 1.9 ^ s ^ 5.2 using tq = 10 nm. These 
values appear realistic given the fact that a 15°C change 
in temperature is enough to change the solubility of CO 2 
by a factor of two. 

It is clear that in order to reliably obtain the depen¬ 
dence of the surface tension on the supersaturation ratio 
one needs to know the contact radius to within a few 
nanometers. In contrast, in order to prove that the solu¬ 
tion is supersaturated and that the surface tension is less 
than the saturated value one does not need to know the 
tip radius precisely. 

On this basis one can conclude that the present analy¬ 
sis of these experimental measurements on nanobubbles 
explicitly confirm what is required by thermodynamics: 
for nanobubble equilibrium the solution must be super¬ 
saturated, and a supersaturated solution has a lower sur¬ 
face tension than a saturated solution— 

In future experimental studies, electron micrography 
or AFM inverse imaging could be used to get an inde¬ 
pendent estimate of the cross-sectional radius of the tip. 
Also, attempts could be made to explicitly control or to 
measure the supersaturation of the solution. 
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